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Abstract 

We used the exact eigenvectors of the generalized Hubbard Hamiltonian solution to 
predict the transition intensities of the well known GR1 and ND1 transitions at the 
neutral and charged vacancies in diamond. In addition to using exact eigenvectors, 
the method of the calculation is more precise than already reported calculations. The 
quantitative results can exlain recent experimental data very good. 
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I. INTRODUCTION 



Coulson and Kearsly introduced molecule concept for vacancy defect in 
diamond. 1 As the results they used the common method for calculation the 
electronic structure of the molecules in quantum chemistry i.e. Molecular or- 
bital and Configuration interaction. 

In essence in this model, molecular orbital method gives correct symmetry and 
spin of the electronic levels of the vacancy and configuration interaction con- 
siders the correlation effects of e-e interaction and gives correct sequence of the 
electronic levels. Molecular orbital method has been the only way for accounting 
the appropriate symmetry and spin of the vacancy since Coulson and Kearsly 
model. 

Already we showed that generalized Hubbard Hamiltonian can also take into 
account correct symmetry and spin of the system in a completely different 
manner. 2 This approach has origin in the solid state physics and pay more 
attention to include symmetry and spin considerations of the system in the 
form of a Hubbard like Hamiltonian. This is in contrast to quantum chemistry 
approaches (molecular orbital theory) which pay more attention to the con- 
struction of the symmetry and spin adapted wavefunctions. 
Both approaches use some parameters (8 parameters for exact evaluation of e-e 
interaction) which can be obtained by theory or semi empirical methods. These 
parameters can be obtained by choosing appropriate atomic orbital function for 
dangling bonds of the vacancy. 

This evaluation is independent of the model of Hamiltonian calculation. These 
parameters are the origin of the one type of approximations which exist in the 
molecular model. These parameters are the only input of the generalized Hub- 
bard Hamiltonian formalism and after evaluating and putting them into the 
computational scheme the exact eigenvalues and eigenstates can be obtained. 
The situation in the Coulson and Kearsly model and its related models 3 which 
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use molecular orbital techniques is different. These models not only use these 
approximate parameters but also apply some approximation in constructing ap- 
propriate CI wavefunctions as the bases for calculating the Hamiltonian. 
Basically the coefficients of the manually constructed CI wavefunctions in these 
models should be calculated by variational principal. Also there should be some 
hypothesis about the most probable configurations in the ground and excited 
electronic state of the system. Additionally there should be an assumption for 
truncation limit in accounting a high number of determinant Slaters in the ex- 
pansion of the CI wavefunctions. 

These assumptions and approximations, results to a approximately CI wave- 
function for the electronic states of the vacancy in Coulson and Kearsly and 
related models. 1 ' 3 In our last communication 2 we also explaineded some ambi- 
guity of the applying such CI method for this problem. 

In summary, it is obvious that the CI wavefunctions in molecular models are 
not the exact eigenvectors of the Hamiltonian of the vacancy system. 
Here we should note to this point that our using bases consist of the atomic 
bases i.e. a, b, c, d in contrast to Coulson and Kearsly and all of the previous 
related models which use molecular symmetric and anti symmetric bases i.e. a, 
t x , t y , t z . These bases are more convenient in explaining physical properties of 
the system. For example ionization of the atoms of the vacancy or possibility 
of the exsistance of the hunds configurations in the ground state can be verified 
very easily. 

Previously we reported 2 the probablity of finding each electronic configuration 
in the exact eigenstates of the generalized Hubbard Hamiltonian for the GR1 
transition in the neutral vacancy. (Fig.l) 

These results and similar results for the ground and excited states of ND1 
transition are free from approximations of the manually constructing CI wave- 
functions. In this method the origin of the error only is in the Hamiltonian 
parameters. 
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FIG. 1: Percentage of each possible electronic configurations in the ground and excited 
states of the neutral vacancy in diamond 



For verifying the validity of these new results we have used them to calculate 
some experimentally observable data. 



II. MODEL 

Similar to Lowther 4 and Mainwood and Stoneham 5 we attempted to calculate 
the relative intensities of the well known GR1 and ND1 transitions in the neutral 
and charged vacancies in diamond. These intensities are calculated directly with 
the resultant eigenvectors of the generalized Hubbard Hamiltonina and can be 
a reliable check for the validity of the physics which is included in the exact 
eigenstates of the model. 

The rate of the electronic dipole transition between an initial state \17j and final 
state is proportional to square of the dipole transition amplitude i.e. 

/(i^/)cx i <^i r i^/>i 2 (i) 

Accordingly for the dipole transition from the ground states of the V° and V~ 
i.e. GR1 and ND1 we have two allowed transitions which their intensities or 
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transition rates are as follow: 

I(GR1) oc | < 4 E | ~r I 1 T > | 2 (2) 

I(ND1) oc | < 4 A | T" | 4 T > | 2 (3) 

The eigenvectors of the generalized Hubbard Hamiltonian have correct degener- 
acy due to symmetry and spin multiplicity of the states. So we obtain two and 
three independent eigenvectors for the ground and excited state of GR1 i.e. l E 
and l T respectively and we will have four and twelve independent eigenvectors 
for the ground and excited states of ND1 i.e. 4 A and 4 T. In the end of this 
paper we have listed ( Table a, b) the expansion coefficients of each of these 
degenerate levels on the S z representation of the solving Fock space bases of 
the Hamiltonian. In the first quantization language this bases are determinant 
Slaters which are correspond to each electronic configuration. For simplicity, in 
this paper we will call these bases as determinant. 

In Table a, b all of the exact eigenstates are orthonormal. An interesting point 
in the Table a is high countribution of the Hunds configurations i.e. (1, 1, 1, 1) 
in the ground state of V° and also the absence of this configuration( model gives 
exactly zero! ) in the the excited state of GR1 . 

At the begining of the calculation we note that this transition amplitudes has 
the same value for the two set of the eigenvectors which transform with a uni- 
tary matrix. So the value of the transition amplitude is the same for S z or S 2 
representation of the Hamiltonian eigenvectors. 

For GR1 transition only the bases or determinant Slaters which have S z — 
enter into the calculation. (Table a) 

36 

\ 1 E,i>=^»)^ (i = l,2) (4) 

3=1 
36 

\ 1 T,i>=Y,^ j (2 = 1,2,3) (5) 

3=1 
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In the above equations i is the degeneracy label of the states and ip/s are the 
S z — bases of the Fock space or determinants. The a* and (3j coefficients are 
according to Table a, b. For GR1 transition we have 6 independent transition 
from the two degenerate l E states to three degenerate l T states. In this work 
we have calculated each transition rate (intensity) separately for each of these 
transitions. 

Here we review in detail the steps of the calculation of the intensity of one 
typical transition e.g. \ l E, 1 > to \ X T,1 >. The amplitude of this transition 
from Eq. (4), Eq. (5) is: 

36 

< x E, 1 | ~r I 1 T, 1 >= ^2 a ^3 <Vi\^\^j> (6) 

In this calculation we have evaluated all of the < <pi | ~r \ tpi > terms up to s 
order, s is the overlap integral of the two adjacent atomic orbital of the vacancy. 

s=<a\b> (7) 

We have calculated the value of s using appropriate Slater type function with 
exponent 1.595 as atomic orbital. The result was 0.16 which agrees with the 
previous models. 1 ' 3 

These types of amplitudes i.e. < ipi \ r \ ^ > are direct expectation values 
of the dipole ~r operator which arise from the same determinant Slaters terms 
of the 1 E and l T eigenvectors in Table a. These four electronic terms can be 
expanded to single particle terms as follow. 

< abed | T" | abed >=< a|7 > |a> + <6|T > |6> + <c|7 > |c> + <(i|7 > |(i> +0(s 2 ) 

(8) 

the terms which are proportional to s 2 are neglected equation in the above 
equation. In addition to these direct terms we have considered the cross terms 
in which the ifi and ipj in Eq. (6) are two different determinant Slaters. The 
main contribution of the cross terms comes from such < ipi \ r \ ipj > terms in 
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which ifi and (fj only differ in one orbital. These terms are proportional to s. 
For example for contribution of the abed and abec determinants we have: 

< abed | ~r \ abec >= s<a\~r\a>+s<b\T > \b>+s<c\~r\c> + <d\~r\c> +0(s 2 ) 

(9) 

In these expansion the sign of the terms should be considered according to the 
permutation of that orbital which is different in two determinants. 
For simplifying the last term in Eq. (9) we have used well known Mullikens's 
approximation: 

g 

<d \ r | c>= -(< d\r\d> + <c\r\c>) (10) 

Now by these rules we can calculate the rate of the transition between each 
degenerate ground and excited states by squaring each transition amplitude. 
For simplifying the writing of the formula we have used the abbreviation: 

~r aa =< a I ~r \ a > (11) 
and therefore the final amplitude of the dipole transition will be in the form: 
< x E, 1 | ~r I 1 T, 1 >= a^aa + a 2 'r'bb + oi^yr cc + a A ~r d d (12) 

The aj are obtained from coefficients which are listed in Table a and also using 
value of s — 0.16. 

Since the ~r aa , ~r bb are not orthogonal and they have equal length, after 
squaring the amplitude to find the transition rates we will have: 

4 4 

K 1 E, 1 | ~r I 1 T, 1 > | 2 = a i + cos(Z109.5°) ^ a^j) \~r aa \ 2 (13) 

i=l i,j=i 

III. RESULTS AND DISCUSSIONS 

The intensities of each six and twelve allowed dipole transitions from 
the two and four fold degenerate ground state of V° and V~ ( GR1 and 
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ND1 ) are summarized in Table I, II. These transitions are calculated by 
the method which was explained in the last section. In the absence of 
L — S coupling, transition between ground and excited states with different 
S z should be zero, therefore as the results of our calculations show, Table 
II should be block diagonal. The total rates of the transitions from each of 
the degenerate ground state of the V° and V~ are summarized in the Table I, II. 





Ground States 


\ 1 E,1> ^E^ > 


Excited states 


S z 








^T.l > 





0.185 


0.051 


| J T,2 > 





0.025 


0.092 


1^,3 > 





0.055 


0.067 


Total Rate: 




0.26 


0.21 



TABLE I: Intensities of the transitions in the unit of |r aa |(|r aa | = f)from the twofold 
degenerate ground state 1 E to threefold degenerate exited state 1 T with S z = for 
the GR1 line. 

As it is seen from the tables the total rates of the transition from each degenerate 
ground state are nearly equal. Physically this sounds good. We expect that in 
an spectroscopic experiment which does not distinguish the degeneracy of the 
ground and excited states, one can not create a change in the relative population 
of the balanced degenerate ground states by illumination of light. 6 
For obtaining the effective transition rate between degenerate ground states we 
should average over the the total transition rates. 6 This arises from this point 
that the probability of the existence of the system in each of the degenerate 
states is equal. In other hand the number of the electrons which populate the 
degenerate ground states are equal. 

Total transition rates of the GR1 and ND1 degenerate ground states are listed 
in the last line of the Table I, II. 
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Ground states \ 4 A, 1 > \ 4 A,2 > | 4 A,3> | 4 A,4> 



Excited states 


Sz 


-3/2 


-1/2 


1/2 


3/2 


| 4 T,1 > 


-3/2 


0.3 











| 4 T, 2 > 


-3/2 


0.31 











| 4 T,3 > 


-3/2 


0.31 











| 4 T,4 > 


-1/2 





0.14 








| 4 T,5 > 


-1/2 





0.24 








i,0 > 


-1/2 


U 


U.oo 


U 


u 


| 4 T,7 > 


1/2 








0.23 





| 4 T,8 > 


1/2 








0.075 





| 4 T,9 > 


1/2 








0.44 





| 4 T, 10 > 


3/2 











0.41 


| 4 T, 11 > 


3/2 











0.32 


| 4 T, 12 > 


3/2 











0.0094 


Total Rate: 




0.92 


0.76 


0.74 


0.74 



TABLE II: Intensities of the transitions in the unit of |r aa |(|r aa | = l)from the fourfold 
degenerate ground state 4 A to twelvefold degenerate exited state 4 T with available 
S z for the ND1 line. 

From Table I we obtain: 

2,3 



| < x E | r I 1 T > | 2 = 1 ]T | <! E,i | t I 1 T,j > | 2 = 0.23|V aa ! 2 (14) 
and from Table II we obtain: 



2 . 



< 4 A\J? \ 4 T > \ 2 = 1-^2 \ < 4 A,i\~r \ 4 T,j > \ 2 = 0.79\^aa\ 2 (15) 



With these values we obtain the relative intensities of the GR1 and ND1 
absorption lines: 

This is in excellent agreement with the recent experimental work on the GR1 
and ND1 transitions intensities. 7 ' 8,9 We have used these values from Davies 10 
who has calibrated them. The experimental data shows that the relative inten- 
sities of the ND1 and GR1 is equal to 3.25 7 ' 8 or 4 7,9 . Our theoretical value is also 
consistant with the experimental error which is reported more recently. 11 The 
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error of our calculation which arises from averaging over the total transition 
rates are 10 percent for ND1 and 11 percent for GR1. These errors should be 
due to approximations in obtaining the Hamiltonian parameters. This means 
that if we were able to estimate exact atomic orbitals of the vacancy and also 
Hamiltonian parameters then we expect that the total rates of the transition be 
equal, (last line of Tables I, II) 

One important question in this problem has been the importance of the over- 
lap effect of the vacancy orbitals. We have recalculated the intensities with 
the neglecting the overlap of the orbitals i.e. s = 0. The results for the GR1 
intensity is 0.23 which is same as Eq. (14) and for ND1 the result is 0.87, 
higher than the value of the Eq. (15). By this assumption we reach to 3.8 for 
the relative transitions. The independent nature of the GR1 transition from 
the overlap parameter s comes from the high contribution of the Hunds state 
in the ground state. As it was shown in Eq. (8) the transition amplitude of 
these type determinants is independent from s. However the value of the s 
depends on the choosing of the atomic orbitals of the vacancy so by changing it 
the Hamiltonian parameters also should be changed. If we were able to include 
semiempirical nature of the generalized Hubbard Hamiltonian parameters to the 
shape of the atomic orbitals then we can be sure that the resultant overlap value 
is more valid. However since the Hamiltonian parameters are from Coulson and 
Kearsly calculation, the value of the s in our calculation are choose to be 0.16. 
Our model is more similar to Mainwood and Stoneham model since their model 
does not need to detail form of atomic orbital i.e. calculation of the | r aa \. Also 
they have attributed the GR1 and ND1 transition only to extra atomic transi- 
tions. They have neglected the overlap of the vacancy orbitals (s = 0) maybe 
due to the order of the error which exist in the approximated wavefunctions of 
the Lannoo model. 3 
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IV. CONCLUSION 



We applied the exact wavefunctions of the generalized Hubbard Hamiltonian 
to the transition intensities of the well known optical absorption lines of the di- 
amond vacancies. This calculation can successfully explain recent experimental 
data. These results clearly show the advantage of using exact wavefunctions of 
the generalized Hubbard Hamiltonian with respect to previous CI wave func- 
tions. For comparing the different wavefunctions of the two model it is useful 
to compare the number of determinants which comes into the calculation of 
the transition rates. For example CI only predicts lxl determinants for ND1 
ground and excited states however here we considered 8x8 determinants for 
the S z = ±| and 48 x 48 for S z = ±| as are reported in Table b. Previous 
molecular orbital approaches were not able to consider different forms of the 
degenerate eigenvectors of the Hamiltonian and this can be another reason fort 
their poor quantitative results. 

The results of the calculation show that the electronic configurations of the 
ground and excited states of GR1 ( Fig. 1) can explain quantitatively experi- 
mental data. This means that the ground state of the neutral vacancy mainly 
comes from the Hunds rule( 60 percent). High contribution of the Hunds states 
in the ground state of the neutral vacancy physically was expected. 
New results of present model show that generalized Hubbard Hamiltonian has 
potential to go toward a quantitative understanding of the point defects problem 
in diamond. 
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Table a. Expansion coeficients of the *E and X T states of the GR1 on 
the S, representation of the Fock space bases (determinants) 
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.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


16 . 


1 /O 
-1/ Z 


aabcc 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


1 7 


-1 /O 

-1/Z 


abccd 


0. 


000 


0. 


289 


0. 


000 





.000 


n nnn 
u . uuu 


0. 


000 


0. 


000 


n 


"3Q 
Z JO 


1 Q 


1 /O 
-1/Z 


aabbc 


0. 


.000 


0. 


000 


0. 


000 





.000 


n nnn 
u . uuu 


0. 


.000 


0. 


.000 


n 

u . 


nnn 
uuu 


19 . 


1 /O 
-1/ Z 


abcdd 





.000 


-0 


.289 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


382 


20 . 


1 /o 
-1/2 


abccd 


0. 


000 


0. 


289 


0. 


000 





.000 


0.000 


0. 


000 


0. 


000 


. 


238 


21 . 


1 /O 
-1/ ^ 


abcdd 





.000 


-0 


.289 


0. 


000 





.000 


0.000 





.000 





.000 


. 


382 


22 . 


1 /O 
-1/ Z 


accdd 





.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


3 
Z D . 


-1 /O 
-1/Z! 


abbcd 


0. 


.000 


-0 


.289 


0. 


000 





.000 


n nnn 
u . uuu 


0. 


.000 





.000 


u 


31 Q 
. JlO 


24 . 


1 /o 
-1/ ^ 


abccd 





.000 


0. 


289 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


238 


25 . 


1 /o 
-1/ Z 


abcdd 


0. 


000 


-0 


.289 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


000 


. 


382 


26 . 


1 /O 
-1/ Z 


ppccg 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 





.000 


. 


000 


07 

z / . 


1 /O 
-1/Z 


DDCQQ 


0. 


000 


0. 


000 


0. 


000 





.000 


n nnn 
u . uuu 


0. 


000 


0. 


000 


u . 


nnn 
uuu 


28 . 


1 /O 
-1/ Z 


bccdd 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 





.000 


0. 


.000 


. 


000 


29 . 


+1/2 


aacaa 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 





.000 


. 


000 


30 . 


+1/2 


aabcc 


0. 


000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


000 


0. 


000 


. 


000 


31 . 


+1/2 


abccd 


0. 


.000 


0. 


000 


-0 


.289 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


"30 
JZ . 


J.1 /o 
+ 1/ z 


aaDDa 


0. 


.000 


0. 


000 


0. 


000 





.000 


n nnn 
u . uuu 


0. 


000 


0. 


.000 


n 
u . 


nnn 
uuu 


33 . 


+1/2 


abcdd 


0. 


.000 


0. 


000 


0. 


289 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


34 . 


+1/2 


aabcd 





.000 


0. 


000 


-0 


.289 





.000 


0.000 





.000 





.000 


. 


000 


35 . 


+1/2 


ajjDcc 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


36 . 


+1/2 


aabcd 


0. 


.000 


0. 


000 


-0 


.289 





.000 


0.000 


0. 


.000 





.000 


. 


000 




J.1 /o 
+ 1/ z 


abbcd 


0. 


000 


0. 


000 


0. 


289 





.000 


n nnn 
u . uuu 


0. 


000 


0. 


000 


n 
u . 


nnn 
uuu 


"3Q 

JO . 


J.1 /o 
+1/Z 


oDDCa 


0. 


.000 


0. 


000 


0. 


289 





.000 


n nnn 
u . uuu 


0. 


.000 


0. 


.000 




nnn 
uuu 


39 . 


+1/2 


,UUJJ 

aDDaa 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


40 . 


+1/2 


abbcd 


0. 


000 


0. 


000 


0. 


289 





.000 


0.000 


0. 


000 


0. 


000 


. 


000 


41 . 


+1/2 


aocca 





.000 


0. 


000 


-0 


.289 





.000 


0.000 





.000 





.000 


. 


000 


42 . 


+1/2 


aapaa 


0. 


000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


000 


. 


000 


A "3 


J.1 /o 
+ 1/ z 


DDCuu 


0. 


.000 


0. 


000 


0. 


289 





.000 


n nnn 
u . uuu 


0. 


.000 





.000 


n 
u . 


nnn 
uuu 


44 . 


+1/2 


abccd 


0. 


.000 


0. 


000 


-0 


.289 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


A C 


+1/2 


abcdd 


0. 


000 


0. 


000 


0. 


289 





.000 


n nnn 
u . uuu 


0. 


000 


0. 


.000 


n 


nnn 
uuu 


46 . 


+1/2 


accdd 


0. 


.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


. 


000 


d 7 


J.1 /o 
+ 1/ z 


aabcd 


0. 


000 


0. 


000 


-0 


.289 





000 


n nnn 
u . uuu 


0. 


000 


0. 


000 


u . 


nnn 
uuu 


48 . 


+1/2 


bbccd 





.000 


0. 


000 


0. 


000 





.000 


0.000 





.000 





.000 


. 


000 


49 . 


+1/2 




0. 


.000 


0. 


000 


0. 


000 





.000 


. 000 


0. 


.000 





.000 


o . 


000 


50. 


+1/2 


bbcdd 


0. 


000 


0. 


000 


0. 


000 





000 


0.000 


0. 


000 


0. 


000 


0. 


000 


51. 


+1/2 


aabbc 





.000 


0. 


000 


0. 


000 





.000 


0.000 


0. 


.000 


0. 


.000 


0. 


000 


52. 


+1/2 


bccdd 


0. 


.000 


0. 


000 


0. 


°°f4 





.000 


0.000 


0. 


.000 


0. 


000 


0. 


000 


53. 


+3/2 


abcdd 


0. 


.000 


0. 


000 


0. 


000 


-0 


.500 


0.000 


0. 


.000 





.000 


0. 


000 


54. 


+3/2 


abbcd 





.000 


0. 


000 


0. 


000 


-0 


.500 


0.000 





000 


0. 


.000 


0. 


000 


55. 


+3/2 


abccd 


0. 


000 


0. 


000 


0. 


000 





.500 


0.000 


0. 


000 


0. 


000 


0. 


000 


56. 


+3/2 


aabcd 


0. 


000 


0. 


000 


0. 


000 


0. 


.500 


0.000 


0. 


000 


0. 


.000 


0. 


000 
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Table b. Expansion coeficients of the A and T states of the ND1 on 



+hft S_ rpnrpspnffll-inn nf t-ho Fnr.lc soarp hasps (Hpl-prmi nani-s) 



S z (Pi | 4 T,5 > | 4 T,6 > | 4 T,7 > | 4 T,8 > | 4 T,9 > | 4 T,10 > | 4 T, 11 > | 4 T,12 > 



1. 


-3/2 


aabcd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


.000 


2. 


-3/2 


abbcd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


3. 


-3/2 


abccd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


000 


4. 


-3/2 


abcdd 


0. 


000 





.000 


0. 


000 


0. 


000 





,000 


0. 


.000 


0. 


.000 


0. 


.000 


5. 


-1/2 


aabcd 


0. 


428 





.193 


0. 


000 


0. 


000 





.000 





.000 





.000 


0. 


.000 


6. 


-1/2 


aabcd 


0. 


428 





.193 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


7. 


-1/2 


aabdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 





.000 


0. 


000 


8. 


-1/2 


aabcd 


0. 


428 





.193 


0. 


000 





000 





.000 


0. 


000 


0. 


000 


0. 


.000 


9. 


-1/2 


aaccd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


10. 


-1/2 


aacdd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


11. 


-1/2 


abbcc 


0, 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


.000 


12. 


-1/2 


aabbd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 





.000 





.000 


0. 


.000 


13. 


-1/2 


abbcd 


-0 


.101 





.373 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


14. 


-1/2 


abbcd 


-0 


.101 





.373 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


15. 


-1/2 


abbdd 


0. 


000 





.000 


0. 


000 


0. 


000 





,000 


0. 


.000 


0. 


.000 


0. 


000 


16. 


-1/2 


aabcc 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


000 


17. 


-1/2 


abccd 


-0 


.259 





.356 


0. 


000 


0. 


000 





.000 


0. 


.000 





.000 


0. 


000 


18. 


-1/2 


aabbc 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


.000 


19. 


-1/2 


abcdd 


0. 


270 





.176 


0. 


000 


0. 


000 





.000 





.000 





.000 


0. 


000 


20. 


-1/2 


abccd 


-0 


.259 





.356 


0. 


000 


0. 


000 





.000 


0. 


.000 





.000 


0. 


000 


21. 


-1/2 


abcdd 


0. 


270 





.176 


0. 


000 





000 





.000 


0. 


000 


0. 


000 


0. 


.000 


22. 


-1/2 


accdd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


23. 


-1/2 


abbcd 


-0 


.101 





.373 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


000 


24. 


-1/2 


abccd 


-0 


.259 





.356 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


.000 


25. 


-1/2 


abcdd 


0. 


270 





.176 


0. 


000 


0. 


000 





.000 


0. 


.000 





.000 


0. 


.000 


26. 


-1/2 


bbccd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


000 


0. 


000 


27. 


-1/2 


bbcdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 





.000 


0. 


000 


28. 


-1/2 


bccdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


.000 


29. 


+1/2 


aacdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


.000 


30. 


+1/2 


aabcc 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


31. 


+1/2 


abccd 


0. 


000 





.000 


0. 


458 


-0 


.131 





.150 


0. 


000 


0. 


000 


0. 


.000 


32. 


+1/2 


aabbd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 





.000 





.000 


0. 


.000 


33. 


+1/2 


abcdd 


0. 


000 


0. 


.000 


0. 


255 


0. 


409 





.135 


0. 


.000 


0. 


000 


0. 


000 


34. 


+1/2 


aabcd 


0. 


000 





.000 


-0 


.239 


0. 


229 





.375 


0. 


.000 





.000 


0. 


000 


35. 


+1/2 


abbcc 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


36. 


+1/2 


aabcd 


0. 


000 





.000 


-0 


.239 


0. 


229 





.375 


0. 


.000 


0. 


.000 


0. 


000 


37. 


+1/2 


abbcd 


0. 


000 





.000 


-0 


.035 


-0 


.311 





.390 


0. 


.000 


0. 


.000 


0. 


000 


38. 


+1/2 


abbcd 


0, 


000 





.000 


-0 


.035 


-0 


.311 





.390 


0. 


000 


0. 


.000 


0. 


000 


39. 


+1/2 


abbdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


000 


40. 


+1/2 


abbcd 


0. 


000 





.000 


-0 


.035 


-0 


.311 





.390 


0. 


.000 





.000 


0. 


000 


41. 


+1/2 


abccd 


0. 


000 





.000 


0. 


458 


-0 


.131 





.150 


0. 


000 


0. 


000 


0. 


.000 


42. 


+1/2 


aabdd 


0. 


000 





.000 


0. 


000 


0. 


000 





,000 


0. 


.000 


0. 


.000 


0. 


000 


43. 


+1/2 


bbcdd 


0. 


000 


0. 


.000 


0. 


255 


0. 


409 





.135 


0. 


000 


0. 


.000 


0. 


000 


44. 


+1/2 


abccd 


0. 


000 





.000 


0. 


458 


-0 


.131 





.150 


0. 


000 


0. 


000 


0. 


.000 


45. 


+1/2 


abcdd 


0. 


000 





.000 


0. 


255 


0. 


409 





.135 





.000 





.000 


0. 


.000 


46. 


+1/2 


accdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


000 


0. 


000 


47. 


+1/2 


aabcd 


0. 


000 





.000 


-0 


.239 


0. 


229 





.375 


0. 


.000 





.000 


0. 


000 


48. 


+1/2 


bbccd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


000 


49. 


+1/2 


aaccd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


000 


0. 


000 


50. 


+1/2 


bbcdd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


.000 


51. 


+1/2 


aabbc 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


000 


0. 


000 


0. 


.000 


52. 


+1/2 


bccdd 


0. 


000 





.000 


0. 


ooj 5 


0. 


000 





.000 





.000 





.000 


0. 


.000 


53. 


+3/2 


abcdd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


.000 


0. 


.000 


0. 


866 


54. 


+3/2 


abbcd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


813 


0. 


074 


-0 


.289 


55. 


+3/2 


abccd 


0. 


000 





.000 


0. 


000 


0. 


000 





.000 


0. 


.471 


-0 


.667 


0. 


.289 


56. 


+3/2 


aabcd 


0. 


000 


0. 


.000 


0. 


000 


0. 


000 





.000 


0. 


.343 


0. 


.741 


0. 


289 
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Table b. (countinued) Expansion coeficients of the A and T states of 



